#include <iostream>
#include <fstream>
using namespace std;

#include "CGravitationalAcretion.h"
#include "../Particle/CParticle.h"  
#include "../SolarNebula/CSolarNebula.h"
 

typedef CParticle<CParticleTraits>	particle_t;
typedef CSolarNebula<particle_t>	solar_nebula_t;
typedef CGravitationalAcretion<solar_nebula_t> gravitational_acretion_t;

int main (int , char * []) {
	cout << "Making nebula" << endl;
	solar_nebula_t  _s ("../SolarNebula/coil100t.txt", 12, 7200, 84, true);
//	solar_nebula_t  _s ("../SolarNebula/iris.txt", 4, 150, 5, true);
//	cout << _s << endl;
	gravitational_acretion_t _g(_s._nebula,_s._size, .01, .1, 10, 20, .5);
	ofstream file("res.txt");
	file << _g << endl;
	file.close();
	system("@start res.txt");

//	for (int i = 0; i < (int)_g._planetesimals.size(); i++)
//		for (int k = 0; k < (int)_g._planetesimals[i]->_container.size(); k++)
//			cout << _g._planetesimals[i]->_container[k] << endl;
	return true;
}